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Abstract 

We present a rational approximation for the Dawson’s integral of 
real argument and show how it can be implemented for accurate and 
rapid computation of the Voigt function at small y << 1. The al¬ 
gorithm based on this approach enables computation with accuracy 
exceeding 10“^^ within the domain 0 < x < 15 and 0 < y < 10“®. 
Due to rapid performance the proposed rational approximation runs 
the algorithm without deceleration. 

Keywords: Dawson’s integral, Voigt function, Faddeeva function, 
complex error function, rational approximation, spectral line broad¬ 
ening 


1 Introduction 

The Dawson’s integral is defined by PElElll] 


2 



0 


*Dept. Earth and Space Science and Engineering, York University, Toronto, Canada, M3J 1P3. 
^Dept. Physics and Astronomy, York University, Toronto, Canada, M3J 1P3. 


1 



where z = x + iy is a. complex argument. This integral is closely related to 
the complex error function, also known as the Faddeeva function [BEIE] 



In particular, comparing equations Q and (|^ we can easily hnd a relation 
between the Dawson’s integral and the complex error function 


F{z) = 




(3) 


As we can see from this identity, the Dawson’s integral is simply a reformu¬ 
lation of the complex error function. 

The Dawson’s integral Q cannot be taken analytically in terms of elemen¬ 
tary functions in a closed form and, therefore, it must be solved numerically. 
Several useful approximations for the Dawson’s integral have been reported 
laEiii]. In this paper we show a new approach that can be successfully 
implemented for efficient computation. 

The complex error function ([^ can be alternatively represented as (see 
Appendix A) 


w{x,y) 



OO 

J exp (—1^/4) exp {—yt) exp (ixt) dt. 
0 


and expressed as a sum |6] 

w{x,y) = K {x,y)+iL{x,y), 2 /> 0, 

where its real and imaginary parts are (see Appendix B) 


and 


K {x,y) 


L {x, y) 




rj dt « 

2/2 + (x - t) 

2/2 + (x - tf 


(4) 


(5) 

( 6 ) 
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respectively. The equation ([^ is known as the Voigt function [SI El El El ED] , 
widely used in many branches of Applied Mathematics [II1E21E3], Physics 
[m IT5l ESI EZl ES ES] and Astronomy [20]. Since the function 0 has no a 
specihc name, it can be regarded as the L-function. 

Due to symmetric properties of the Voigt function 

K (x, - \y\) = K (-x, - \y\) = -K {x, \y \), 

it is sufficient to consider only the P* and quadrants in order to cover 
the entire complex plane. Therefore, we will imply further that y > 0. 

In our recent publication [10] we have shown that a new sampling method¬ 
ology based on incomplete expansion of the sine function IZH leads to a ratio¬ 
nal approximation of the Voigt function for rapid and accurate computation 
(see also the Matlab source code in m) 

’ ^ + 2/^m (2/^ - + (a;2 + y'^f (7) 

^ K {x,y) ^ K{x,y + ‘;/2 ), 


where 


^ 0t(m- 1/2) A ,2/4_„2;,2 / 7r(m-l/2) (nh + <?/2) 


Pm 


/vr (m — 1/2) 
V 2mmax^ 


2 


5 


1m 



N 

V cos 


n=-N 



1/2) {nh + <?/2) 

^maxh 


h = 0.293, TTi max = 12, ^ = 2.75 and N = 23 (for higher accuracy the 
coefficients h and mmax can be taken as 0.25 and 16, respectively). This 
rational approximation is faster than the Weideman’s rational approximation 
[22] by factor greater than two (see ra for details). Furthermore, with the 
same number of summation terms it is by several orders of the magnitude 
more accurate than the Weideman’s rational approximation in the domain 
of practical interest 0 < x < 40, 000 and 10“ ^ <y < 10“2 [1^123] required 
for applications using the HITRAN molecular spectroscopic database [2T| . 
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In general, the series approximation ([^ provides accurate results while y > 
10“®. However, its accuracy deteriorates with decreasing parameter y and, 
consequently, it cannot cover a narrow band region y < 10“® along x-axis. 
It should be noted that the accuracy deterioration of the Voigt function at 
small ?/ << 1 is a very common problem in most known approximations (see 
for example [3 EH 12S] ) ■ 

Although the computation of the Voigt function ([^ at ?/ < 10“® is re¬ 
quired relatively rare in practice, it, nevertheless, has to be taken into ac¬ 
count in algorithmic implementation. In this work we propose a rational 
approximation for the Dawson’s integral of real argument and show how its 
implementation can resolve effectively such a problem in computation of the 
Voigt function that occurs at small y « 1. 


2 Results and discussion 


2.1 Derivation 


As we have shown recently, the exponential function can be expanded as a 
series [10] 


-I '^max 

exp (-x^) ~ 2 £ 

m=l 


+ (x -|- Bra ^ + {—X -|- Bm 

- Cl, - {x + iql2f C^-(-x + k/2)^ . 


( 8 ) 


where 


■^m O'm 


y/n{m- 1/2) 
21711^ h 


N 

E 

n=-N 


/4_ 

sm 


7T {m — 1/2) {nh -|- <^/2) 


rur 


rh 


Bm — fdm 



N 


cos 

n=-N 



1/2) {nh + c;/2) 

'^raa.y.h 


f 7r(m- l/2) y 

\ 2?7lniax^ / 


This equation has been used to obtain the rational approximation ([^ ra- 
We may also utilize it for the L-function by substituting the exponential 
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function approximation ([^ into equation ([^. This leads to the following 
integral 


L{x,y) « 


OO 

I p I ^^max 

/ X — t 

27r y y^ + {x- tf ^ 

— OO 


+ (i + k/ 2) -Bjn 


dm + i~t + *‘?/2) Bm 


(9) 


The integrand of this integral is analytic everywhere over the entire complex 
plane except 2 + drumax isolated points 


{x iy^ X + iy, Cm C*m C'm T Cm T f?/2} ; 

771 G {1, 2, 3, ... Ulmax} 


where we can observe the singularities. In order to integrate equation (9), we 
may choose conveniently a contour Cccw in counterclockwise (CCW) direction 
on the upper-half complex plane as a semicircle centered at the origin with 
inhnite radius. Consequently, since the domain enclosed by contour Cccw 
covers only the half of complex plane, the number of the isolated points is 
reduced twice 


tr {x -|- iy^ Cm T ^* 5 / 2 , Cm T 5 ni G {1, 2, 3, ... Ulmax 
Using the Residue Theorem’s formula now 


}• 


H-2mn 


— (j) f{t)dt= Res [/ (t ), U], 


Cc. 


r=l 


where / (f) is the integrand of integral ([^, we obtain a rational approximation 
for the L-function (1^ as follows 


^ix,y) = Y 


X [ 2 amy + 7m (2:^ + y‘^ - / 3 m)] 


13m + 2/3m {y‘^ - x"^) + (a;2 + 7/2)2 
^ L{x,y) ^ X{x,y + ^/2). 

Consequently, the complex error function can be calculated as 
w (x, y) ^ K (x,y + q/2) + iX (x, y + (,/2). 


( 10 ) 


( 11 ) 


The computational testing we performed with equation 01 shows that its 
imaginary part X{x,y + <;/2) remains always accurate regardless how small 
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the parameter y is taken. Therefore the computation of the Voigt function at 
small y « 1 is our main objective. We may attempt to resolve this problem 
by using the following approximation that has been reported in our recent 
work |26] 


w (x, y « 1) 


1 + 


2F (x) 


_ ^2ixy 


X 


where the Dawson’s integral of real argument is denoted as F [x) = F {x,y = 0). 
Taking the real part of the equation above yields the approximation of the 
Voigt function 

2 _ 2 

K {x,y « 1) Rs cos {2xy) -^ [y sine {2xy) — F (x) sin {2xy)] , (12) 

where we imply that the sine function satishes 

{sine {2xy 7 ^ 0) = sin {2xy) / {2xy ), sine {2xy = 0) = 1} . 

According to equation ([^ the Dawson’s integral of complex argument can 
be written now as 

F (x, y) ~ ^ - [k (x, y + ?/2) + i\ (x, y + q/2)\ | 

and since at K {x,y = Q) = exp {—F) (see Appendix B), from equation above 
we get an approximation for the Dawson’s integral of real argument 


TT . 


F{x,y = Gi)p^ —\{x,q/2). 


(13) 


Lastly, using the series approximation (10) the equation (13) can be expressed 
as given by 


TTlmax 

F{x)^^^ 


X [am(. + 7m + <tV4 - /^m)] 


(14) 


/^m + (?V4 - x'^) + (a;2 + <^2/4)^ 

Figure 1 shows the difference £ (f) between the original Dawson’s integral 
of real argument and its rational approximation (14). As we can see from this 
hgure, despite only 12 summation terms involved in the series approximation 
(14), the difference e{t) remains within a narrow range ±7 x 10“®. This 
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Fig. 1. The difference e (t) at <^ = 2.75, N = 23, mmax = 12 and 
h = 0.293. 


confirms a rapid convergence of equation (14) that makes it suitable for 
practical applications. Therefore, we can use it effectively as a supplement 
for computation of the Voigt function at small y << 1. 

2 

Since at ?/ << 1 the following functions can be approximated as ~ 1, 
cos ( 2 x 1 /) ~ 1, sine (2X//) ~ 1 and sin (2x//) 
signihcantly simplified as follows 


2xy, the equation (12) is 


K (x, y « 1) e" 


—^ [1 — 2xF (x)] 

VTT 


and according to equation (13) it can also be represented in form 


K (x, y « 1) e ^ [l — \/ttx\ (x, <^/2)] . 


(15) 


The computational time of the approximation (15) is mostly taken by 
A-function. Since the approximations ([^ and (10) are about equally rapid, 
the approximation (15) is also as fast as the approximation ([^. Thus, 


we can 


see that implementation of the Dawson’s function of the real argument (14) 


does not decelerate the computation of the Voigt function at small y «1. 
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2.2 Implementation 


When the input parameters x and y are large enough, say if the condition 
\x + iy\ > 15 is satished, many rational approximations become effective 
for computation. For example, the Gauss-Hermite quadrature, the Taylor 
expansion series jHI or the Laplace continued fraction ISIl [2S1 ESI can be 
used. Therefore, despite that the rational approximations Q and (10) can 
cover accurately the entire domain of practical interest 0 < x < 40, 000 and 
10-^ < 2/ < 102 [ig |23] required for applications using the HITRAN molec¬ 
ular spectral database [21], we may restrict them only within the domain 
|a; -h iy\ < 15 that is considered the most difficult for rapid and accurate 
computation. Thus, the computation of the Voigt function inside the do¬ 
main \x + iy\ < 15 can be implemented in accordance with the following 
scheme 


K{x,y + <;/2), 


\x + iy\ < 15n 10"® < y < 15 


^ ^ e ^ [l — (x, ?/2)] , |x -|- iy| < 15 n 0 < y < 10 


2.3 Error analysis 

Let us dehne the relative error as 


A = 


K (x,y) - Kref. {x,y) 


Kref. {X,y) 


where Kref. {x, y) is the reference, to quantify the accuracy of the approx¬ 
imation (15). The highly accurate reference values can be generated, for 


example, by using the Algorithm 680 in, the Algorithm 916 or a new 
algorithm descried in the recent work [32] . 

Figure 2 show the logarithm of the relative error log^gA for the Voigt 
function inside the domain 0 < x < 15 and 0 < y < 10“®. As we can see from 
this hgure, despite that only 12 summation terms involved in computation. 


the accuracy of approximation (15) in this domain is better than 10 


1-10 


Thus we can see that the approach based on implementation of the ratio¬ 


nal approximation for the Dawson’s integral of real argument (14) is efficient 
and practical. 









1.x 



Fig. 2. Logarithm of the relative error log^gA for the Voigt function 
approximation (15) inside the domain 0<a;<15n0<?/< 10“® 
computed at c = 2.75, N = 23, rrtmax = 12 and h = 0.293. 


3 Conclusion 


A rational approximation (14) for the Dawson’s integral of real argument is 
presented. It can be utilized as an efficient supplement for accurate com¬ 
putation of the Voigt function at small y « 1. In particular, this ap¬ 
proach enables computation with accuracy exceeding 10“^*^ inside the domain 
0 < X < 15 and 0 < ?/ < 10“®. Due to rapid performance the algorithmic 
implementation of the rational approximation (14) does not decelerate the 
computation of the Voigt function. 
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Appendix A 

There is a simple proof of equation (|^. The equation (|^ can be rearranged 
in form 


wiz) = e H- iz 


TT 


and since 


e f exp (—f^/4) cos (zt) dt, 


TT 


0 

OO 


F (z) = - / exp (—f^/4) sin (zf) dt, 


we can write 


w 


{z) = — 1 = / exp (—1^/4) (cos (zt) + isin ( 2 ;^)) dt 


= / exp (—f^/4) exp {izt)dt 


or 




(x, y) = / exp (—f^/4) exp (i (x + iy) t) dt 


= / exp (—1^/4) exp (—|/t) exp [ixt) dt. 


This completes the proof. 


Appendix B 

Due to discontinuity at t = x the functions (@ and @ cannot be dehned 
rigorously aty = 0. However, since the following limits exist (see for example 
[3 El El] for details) 

2 

lim K (x, y) = e~^ , 


10 


and 


2 


F{x,y = 0) 


lim L (x, y) 

y^o ^ Wti 


we may imply that K {x,y = 0) and L{x,y = 0) are defined and eqnal to 
e“* and -^F {x, y = 0), respectively. 
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